749 

UNCLASSIFIED 


INFINITE  ELEHENTS  FOR  THO-DINENSIOWH.  FLUID-STRUCTURE 
INTERACTION  PROBLEHS(U>  DAVID  H  TAVLOR  NAVAL  SHIP 
RESEARCH  AND  DEVELOPNENT  CENTER  BET. .  E  A  SCHROEDER 
FEB  BE  DTRERDC-BE/aai  F/G  13/19. 1 


1/1 


L'NCLASSIK'  ED 

(  I  ASSif iCATiON  DE'viiS  PAut 


PtPO^’'  SECijRlTY  ClASSTICATION 
I  .'.'CLASSIFIED 


/  '  /  L  ,  /  /  . 


REPORT  DOCUMENTATION  PAGE 


’b  RESTRICTIVE  MARKINGS 


.',1  SftEH.TV  Cl  ASS-fiCATION  AUTHORITY 


.p  .'■(  AS^  ‘  CAUON  DOWNGRADING  SCHEOUIE 


;  .  •  .'R'.- ORoAN  i'AT  iTN  REPORT  NUMaER(S) 

:  NSR!)E:-Sb/0O  L 


)  DISTRIBUTION 'AVAILABILITY  OE  REPORT 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION 
IS  UNI.lMll'EI). 


S  MON' TOR. NO  ORGANISATION  REPORT  NUM8£R(S) 


•  NA'.'i  'I  PI  RTORViNG  organization  6b  OFfICE  SYMBOl  la  NAME  OF  MONiTORiNG  ORGANIZATION 
A.  lavlor  Naval  Ship  {if  applicable^ 

.s'.  Contor  Code  1844 


^  v..fv  SfJfP  /Ip  CiX^e) 

,..si!a,  Maryland  d()084-5000 


'b  ADURf  SS  (C/fv  SfAfe,  aiid  /IPCode) 


UNO  sPONSOR'Nl 


8b  office  symbol  9  PROCURFMEN'  'NSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 


>0  SOURCE  OF  FUNDING  NUMBERS 


61153N 


PROlECT 

TASK 

NO 

NO 

15321 

WORK  UNIT 
ACCESSION  NO 


Svcu'ify  Odisificafion) 

INI  IF.  ELEMENTS  FOR  TWO-DIMENSIONAL  FLUID-STRUCTURE  INTERACTION  PROBLEMS 


■..T.  A  U'ORlS) 

LrotJer,  Frwin  A. 


'  !b  Time  covered 
FRCTM  TO 


14  DATE  Of  REPORT  {Year.  Month.  Day)  IS  PAGE  COUNT 

1986,  February  19 


’8  SLiBiECT  [RMS  (Continoe  on  reverse  if  necessary  and  identify  by  bloclc  number) 
Infinite  Elements 
Fluid-Structure  Interaction 


'  Confinue  on  referee  if  necessary  and  identify  by  blexfc  number) 

Ihis  report  de.scribes  an  effective  method  for  using  two-dimensional  infinite 
dements  to  compute  acoustic  or  magnetic  fields  in  the  unbounded  fluid  region 
surrounding  a  submerged  vehicle.  In  this  method,  finite  elements  represent  the 
Is  uiided  rej’.ion  containing  the  vehicle  and  may  also  be  used  to  represent  a  layer  of 
ilaid  surrounding  the  vehicle.  Infinite  elements  are  used  to  represent  the  unbounded 
oterior  region.  Since  infinite  elements  are  not  bounded,  their  shape  functions  are 
,  hosr-n  to  contain  decay  factors  to  produce  convergent  integrals.  If,  from  physical  or 
I  Lher  considerations,  the  order  of  decay  of  the  solution  as  the  radius  increases  is 
known,  ini inite  elements  should  be  chosen  with  the  same  order  of  decay.  The  results 
ouiained  in  this  study  were  found  to  be  within  2^^rcent^when  the  decay  factor  of 
infinite  elements  matched  that  of  the  solution.  However,  for  other  problems,  the 
oidiT  of  decay  of  the  solution  may  not  be  known  In  advance,  and,  therefore,  it  may  not  — ^ 

(Continued  on  reverse  sideic;^' 


,  ■  Cl’l  .YV  AIL  AB  I  V  DF  ARSTRAC' 

/.SLi'F  Ml  UNL  M'Tf  D  LI]  SAME  AS  »PT  Q  OTIC  USERS 


,  M.  SVfiNSiHl  f  .NDiV'DiJAl 

Lrwin  A.  SclircDeder  _ 


DDfORM  1473  AiTVAR  8!  APR  I'biton  may  be  wrI  vntl  Ounausted 

All  Other  editions  jre  obsolete 


21  abstract  SECURITY  CLASSIf (CATION 

□  one  USERS  UNCLASSIFIED  _ 


ZZb  TELEPHONE  (Inc/uc/e  Area  Code)  ZZc  OFFICE  SY.vIBOL 

(202)  227-1660  Code  1844 


1  may  be  vv'd  vnt'l  o«nausted  SECURITY  Cl  ASSIF  iC  A  TiON  OI  TmS  PA(,E 

er  edit.ony  jie  oDsiilete  UNCLASSIFIED 


_ UNCLASSIFIED _ 

iFCURITY  CLASSIFICATION  OF  THIS  PAGE 


(Block  19  continued) 

be  possible  to  match  the  two  rates  of  decay.  For  such  cases,  the  errors  were 
found  to  be  as  large  as  20  percent.  In  such  situations,  a  layer  of  finite  fluid 
elements,  two  elements  thick  around  the  structure,  reduced  the  errors  to  less 
tlian  J.5  percent  for  the  modes  and  decay  factors  tested. 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


TABLE  OF  CONTENTS 


LIST  OF  FIGURES . 

ABSTRACT . 

ADMINISTRATIVE  INFORMATION . . 

I NTRODUCTION . 

NATURAL  FREQUENCIES  OF  A  SUBMERGED  CYLINDER . . 

KINITF.  AND  INFINITE  ELEMENT  MODEL  OF  A  FLUID-STRUCTURE  SYSTEM 

INFINITE  ELEMENT  FORMULATION . 

RESULTS  OF  THE  TEST  PROBLEMS . 

i'ISCUSSION  OF  RESULTS . . 

REFERENCES . 


LIST  OF  FIGURES 


!  An  Infinite  Element  and  its  Parent  Strip . 

Numinique  Mapping  for  Infinite  Elements . . 

I  --  Use  of  Finite  Elements  to  Avoid  Nonunique  Mapping . 

.  •  Gonparison  of  Relative  Errors  for  16  and  32  Sectors  per 

Ouadrant,  Using  Two  Rings  of  Fluid  Finite  Elements . 

1  -  Comiiari  son  of  Relative  Errors  for  Infinite  Element, 

and  Finite  and  Infinite  Element  Solutions . . . 

Accession  For 

NTIS  GFA&I 
DTIC 

Unnunoum  ''d 

By _ _ _ 

Distribi't,  I  ■■/ 
Aval  lab  ii  ■  * -,■ 


ABSTRACT 


This  report  describes  an  effective  method  for  using 
two-dimensional  infinite  elements  to  compute  acoustic  or 
magnetic  fields  in  the  unbounded  fluid  region  surrounding 
a  submerged  vehicle.  In  this  method,  finite  elements 
represent  the  bounded  region  containing  the  vehicle  and 
may  also  be  used  to  represent  a  layer  of  fluid  surrounding 
the  vehicle.  Infinite  elements  are  used  to  represent  the 
unbounded  exterior  region.  Since  Infinite  elements  are 
not  bounded,  their  shape  functions  are  chosen  to  contain 
decay  factors  to  produce  convergent  Integrals.  If,  from 
physical  or  other  considerations,  the  order  of  decay  of 
the  solution  as  the  radius  Increases  is  known,  infinite 
elements  should  be  chosen  with  the  same  order  of  decay. 

The  results  obtained  in  this  study  were  found  to  be  within 
2  percent  when  the  decay  factor  of  infinite  elements 
matched  that  of  the  solution.  However,  for  other  problems, 
the  order  of  decay  of  the  solution  may  not  be  known  in 
advance,  and,  therefore,  it  may  not  be  possible  to  match 
the  two  rates  of  decay.  For  such  cases,  the  errors  were 
found  to  be  as  large  as  20  percent.  In  such  situations, 
a  layer  of  finite  fluid  elements,  two  elements  thick  around 
the  structure,  reduced  the  errors  to  less  than  3.5  percent 
for  the  modes  and  decay  factors  tested. 


ADMINISTRATIVE  INFORMATION 

This  work  was  performed  at  the  David  W.  Taylor  Naval  Ship  Research  and 
Development  Center  (DTNSRDC),  and  was  funded  by  the  Naval  Sea  Systems  Command 
(S05R24)  under  Task  Area  SROIAO^OI,  Task  15321,  and  Work  Unit  18A0-040. 

INTRODUCTION 

To  compvite  magnetic  or  acoustic  fields  about  submerged  vehicles,  the  vehicle, 
surrounding  water,  and  interactions  between  them  must  be  modeled.  The  field  and 
the  structure  can  each  be  represented  by  finite  elements  and  the  representations 
coupled  using  a  method  developed  by  Zienkiewlcz  and  Newton. If  the  vehicle  is 
submerged  in  the  sea,  the  extent  of  the  water  surrounding  it  is  so  great  that  the 
region  containing  the  water  is  best  represented  as  having  infinite  extent.  If 
finite  elements  are  used  to  solve  such  a  problem,  one  can  model  the  structure  and 
a  reasonable  extent  of  the  finite  region  of  the  medium  around  the  structure.  This 
leaves  the  problem  of  accounting  for  the  infinite  region  surrounding  the  finite 


*A  complete  listing  of  references  is  given  on  page  15. 


region  modeled  with  finite  elements.  Various  methods  have  been  proposed  for  repre¬ 
senting  the  surrounding  Infinite  region.  Several  examples  are: 

•  The  surrounding  region  can  he  truncated  by  modeling  a  large  part 
of  It  with  finite  elements  and  applying  boundary  conditions  that 
approximate  the  effects  of  the  remainder  of  the  infinite  region. 

•  The  solution  can  be  expanded  as  a  series  of  analytic  functions  in 
the  surrounding  region  with  the  coefficients  Introduced  as  unknowns. 

•  The  effects  of  the  surrounding  region  can  be  represented  by  an 
Integral  equation  on  the  boundary  of  the  region  surrounding  the 
finite  elements. 

•  The  surrounding  region  can  be  modeled  with  Infinite  elements. 

Zienkiewlcz  et  al.2  compare  these  methods  and  find  that  each  has  certain  advantages 
and  disadvantages. 

In  this  report,  infinite  elements  are  used  to  represent  the  surrounding  infi¬ 
nite  region.  An  Infinite  element  represents  a  sector  of  the  infinite  region  extend¬ 
ing  from  the  boundary  of  the  finite  region.  In  both  finite  and  infinite  elements, 

Che  unknown  function  Is  approximated  by  shape  functions  and  a  functional  involving 
the  shape  functions  Is  integrated  over  the  area  or  volume  of  the  element.  To  obtain 
convergent  Integrals  for  infinite  elements,  one  of  two  schemes  is  used.^  The  first 
incorporates  decay  factors  In  the  shape  functions  that  vary  in  the  direction  that 
extends  to  Infinity.  The  other  scheme  maps  the  infinite  element  into  a  standard 
s,('!ire  or  cube.  This  mapping  also  requires  decay  factors  that  compress  the  Infinite 
('lement  Into  a  finite  region.  The  decay  factors  may  take  one  of  several  forms;  they 
may  decrease  exponentially  with  the  distance  r  from  a  fixed  point,  or  decrease  as  a 
power  of  r.  If  a  solution  has  a  component  that  decreases  as  1/r  as  r  tends  to  infin¬ 
ity,  usually  that  component  will  dominate  for  very  large  values  of  r.  However,  com¬ 
ponents  that  are  significant  for  moderately  large  r  effect  the  accuracy  obtained 
with  infinite  elements.  For  some  problems,  the  rate  of  decay  of  the  predominant 
component  of  the  solution  as  r  tends  to  Infinity,  may  be  deduced  a  priori  from 
physic, il  or  other  considerations.  For  such  problems,  one  could  choose  a  decay 
factor  that  matches  the  behavior  of  this  component  as  r  approaches  infinity.  It 
would  he  expected  that  this  choice  would  produce  the  most  accurate  results.  How- 
evir,  the  solution  may  have  more  significant  components  than  the  infinite  element 
tias  decay  factors,  or  It  may  have  significant  components  that  are  not  known 


a  priori.  That  Is,  the  solution  may  have  significant  components  that  df'creasi*  as 
l/r®  for  one  or  more  values  of  m,  where  the  factors  1/r™  are  not  represented  by 
the  decay  factors  of  an  infinite  element.  Thus,  the  question  arises  of  how  well 
an  Infinite  element  will  approximate  a  solution  that  has  components  that  decrease 
with  orders  different  from  those  represented  by  the  decay  factors. 

This  report  gives  results  for  the  two-dimensional  problem  of  an  Infinitely 
long  cylindrical  shell  submerged  in  an  infinite  acoustic  fluid.  The  ultimate 
interest  is  in  three-dimensional  problems,  but  a  computational  solution  using 
infinite  elements  for  the  two-dimensional  problem  is  easier  to  develop  and  an 
exact  solution  is  available  for  comparison.  A  two-dimensional  problem  is  quite 
different  in  character  from  a  corresponding  three-dimensional  problem.  Rather 
than  simulating  the  complete  three-dimensional  solution  from  the  two-dimensional 
solution,  this  study  examines  the  performance  characteristics  of  two-dimensional 
infinite  elements  and,  based  on  this,  conjectures  the  performance  of  three- 
dimensional  infinite  elements.  This  is  done  for  two  configurations  and  several 
formulations  of  the  infinite  elements.  The  reported  results  show  that  the  approx¬ 
imate  solutions,  produced  when  only  infinite  elements  are  used  to  model  the  fluid, 
liave  small  errors  if  the  decay  factors  match  the  known  rate  of  decay  of  the  exact 
solution.  Also,  the  results  show  that,  for  cases  when  decay  factors  for  the  solu¬ 
tion  and  tlie  infinite  elements  did  not  match,  the  computed  solutions  were  found  to 
be  quite  accurate  if  a  layer  of  the  fluid  surrounding  the  vehicle  was  modeled  by 
finite  elements. 

The  advantages  of  using  infinite  elements  are;  (1)  they  maintain  symmetry 
and  reasonable  bandwidths  in  the  matrices  produced,  and  (2)  it  is  easy  to  implement 
the  method  using  available  finite  element  programs.  On  the  other  hand,  some  formu¬ 
lations  of  infinite  elements  Introduce  more  degrees  of  freedom  than  do  competing 
metnods.  In  some  cases  extra  effort  must  be  made  to  ensure  element  shapes  that 
provide  unique  mappings.  So,  for  some  problems,  these  characteristics  may  require 
additional  effort  in  producing  the  numerical  model. 


NATURAL  FREQUENCIES  OF  A  SUBMERGED  CYLINDER 

We  considered  the  problem  of  an  infinitely  long  cylindrical  shell  submerged 
in  an  acoustic  fluid.  The  problem  is  to  determine  the  natural  frequencies  of  the 
system  consisting  of  the  vibrating  structure  interacting  with  the  surrounding 
f  luid. 

The  axis  of  the  cylindrical  shell  is  aligned  with  the  z-axls.  The  ni  tion  of 
the  cylinder's  surface  and  the  acoustic  pressure  are  assumed  to  be  independent  of 
the  z-coordinate .  Since  it  is  assumed  that  there  is  no  variation  of  shell  dis¬ 
placements  or  fluid  pressure  in  the  z-directlon,  a  two-dimensional  problem  can  be 
obtained  by  projecting  the  cylinder  and  the  region  containing  fluid  into  the 
X, y-plane. 

FINITE  AND  INFINITE  ELEMENT  MODEL  OF  A  FLUID-STRUCTURE  SYSTEM 

Natural  frequencies  of  the  fluid-structure  system  are  computed  using  a  NASTRAN 
model  containing  both  finite  and  infinite  elements.  Equations  associated  with  the 
cylindrical  shell,  the  acoustic  fluid,  and  the  coupling  between  them  are  combined 
to  obtain  a  matrix  equation  for  the  coupled  system.  For  an  assumed  time  harmonic 
solution,  the  matrix  equation  becomes  the  equation  for  an  eigenvalue  problem  for 
determining  the  natural  frequencies  of  the  coupled  system.  This  work,  extends  and 
modifies  the  work  of  Schroeder  and  Marcus^  by  Introducing  Infinite  elements  and 
using  the  added  mass  approximation. 

In  the  two-dimensional  finite  element  model,  the  cylindrical  shell  is  repre¬ 
sented  by  a  ring  of  one-dimensional  bar  elements.  The  fluid  region  is  the  region 
of  th(’  x,y-plarie  outside  the  ring,  and  is  modeled  using  two-dimensional  membrane 
finite  elements  by  an  analogy  for  the  Laplace  operator  and  two-dimensional  infinite 
e lements. 

Standard  finite  element  modeling  procedures  produce  the  structure  matrix 
equation  associated  with  the  ring^ 


Mil  +  Ku  =  f 


where  M  =  mas.s  matrix 

K  =  stiffness  matrix 

u  =  displacement  vector 

f  =  force  vector  acting  on  the  ring 


4 


The  method  of  analogies  uses  structural  membrane  elements  to  produce  a  matrix 
equation  for  the  acoustic  fluid.  The  method  of  analogies^  is  applied  by  giving 
Young's  modulus  E  the  value  unity,  and  the  shear  modulus  G  the  value  10^.  The 
differential  equation  for  the  acoustic  pressure  P  is  the  wave  equation 

-  c2v2p  =  0 

where  c  is  the  speed  of  sound  in  the  fluid. 

If  the  acoustic  pressure  is  harmonic  in  time,  that  is,  P(x,t)  =  p(x)e^‘'^t 
where  x  is  a  point  in  the  x,y-plane,  the  wave  equation  becomes 

k.2p  +  72p  =  Q 

where  k  =  iV c .  Using  the  method  of  analogies  with  the  finite  element  approximation 
produces  the  acoustic  pressure  matrix  equation 

k2Qp  +  Hp  =  0 

The  structure  and  acoustic  pressure  matrix  equations  are  coupled  by  inter¬ 
actions  at  the  fluid-structure  Interface.  The  acoustic  pressure  acting  on  the 
Interface  produces  a  force  vector  (with  the  convention  that  forces  acting  outward 
from  the  surface  are  positive)  so  the  matrix  equation  for  the  structure  becomes 

Mil  +  Ku  =  -Ap 

where  A  is  a  matrix  whose  entries  reflect  the  area  of  elements  on  the  interface . 

See  Schroeder  and  Marcus'*  for  details  of  the  analysis  of  the  fluid-structure 
coupling.  The  matrix  A  may  be  given  by  either  a  consistent  or  a  lumped  formula¬ 
tion.  Since  the  boundary  in  this  problem  is  an  arc  of  a  circle,  the  consistent 
form  is  used  as  it  is  easy  to  compute.  The  entries  in  the  area  matrix  are 

Ai.i  =  fniNiFjdS 

S 


where  =  the  shape  function  associated  with  the  i-th  structural  degree  of 
freedom 

S  =  the  fluid-structure  interface 

nj^  =  the  component  of  the  unit  normal  to  S  in  the  direction  of  Che  i-th 
translational  degree  of  freedom 

Fj  =  the  shape  function  associated  with  the  j-th  fluid  degree  of  freedom 
The  motion  of  the  structure  also  affects  the  acoustic  pressure  field  in  the 
fluid.  The  acceleration  iin  of  the  surface  of  the  structure  normal  to  the  inter¬ 
face  between  the  fluid  and  the  structure  is  related  to  the  gradient  of  the 
acoustic  pressure  field  by  the  equation 

3p/3n  =  -PUn 

where  P  is  the  density  of  Che  fluid.  This  relation  is  added  as  a  boundary  condi¬ 
tion  and  produces  a  force-type  terra  in  the  equation  for  the  acoustic  press  .re 
field.  Passing  to  the  finite  element  matrix  equation  results  in  the  addition  of 
an  area  matrix  term 

k^Qp  +  Hp  =  -pA'^ii 


The  matrix  A  is  the  same  area  matrix  that  appears  in  the  matrix  equation  for  the 
structure. 

For  low  frequencies,  the  added  mass  approximation  can  be  used.  This  approxi- 
:n,ition  accounts  for  only  the  inertial  effects  of  the  surrounding  fluid.  If  the 
frequency  oj  is  small,  then  k  is  small  and  the  first  terra  in  the  preceding  equation 
is  negligible. 

Combining  the  matrix  equations  for  the  structure  and  the  acoustic  field 
produces  Che  system  of  equations 


Mli  +  Ku  +  Ap  =  0 


-PA'^u  +  Hp  =  0 


The  system  of  two  equations  is  reduced  to  one  equation  by  eliminating  p 


Mil  +  Ku  +  pAH~'a'^u  =  0 


In  this  equation,  the  term  PAH“*A^  is  a  coefficient  of  the  second  time  derivative 
of  the  displacement  of  the  structure,  and  thus  acts  like  a  mass  added  to  the  mass 
of  the  structure  and  is  called  the  "added  mass." 

Assuming  a  solution  of  the  form  u(x,t)  =  results  in  a  matrix  eigen¬ 

value  problem  tliat  is  solved  for  the  natural  frequencies  of  the  fluid-structure 
s  V  s  t  e  m 

dec[(ij2(M  +  pAirUT)  ^  k]  =  0 

INFINITE  ELEMENT  FORMULATION 

The  infinite  elements  used  here  are  of  the  decay  factor  type^;  their  formula¬ 
tion  is  an  extension  of  the  formulation  of  finite  elements.  Each  infinite  element 
represents  a  sector  of  the  infinite  region  radiating  outward  from  the  boundary  of 
tile  region  represented  by  the  finite  elements.  The  difference  between  an  infinite 
element  of  this  type  and  a  finite  element  is  that  the  infinite  element  represents 
an  infinite  region  and  the  shape  functions  associated  with  the  infinite  element 
incorporate  decay  factors,  factors  that  decrease  to  zero  fast  enough  to  produce 
convergent  Integrals  over  the  infinite  region. 

The  formulation  of  these  infinite  elements  is  similar  to  that  of  isoparametric 
finite  elements.  The  region  represented  by  one  infinite  element  is  mapped  onto  a 
standard  infinite  strip.  Figure  1  shows  the  mapping  and  the  standard  strip.  The 
functions  that  map  the  infinite  element  into  the  standard  strip  are 

X  =  (x[  +  ^n|)Lj(ii)  +  (x2  +  Cn2)L2(h) 

where  n(  are  unit  vectors  normal  to  the  boundary  between  the  finite  and  infinite 
elements  at  the  points  xj  (i  =  1,2).  The  pressure  field  is  approximated  by 


P  =  PiH(^)Li(n)  +  p2H(^)L2(h) 


Figure  I  -  An  Infinite  Element  (a)  and  its  Parent  Strip  (b) 


where 


Li(n)  =  (1  -  n)/2,  L2(n)  =  (1  +  n)/2 


and 


for  R  =  R(n),  the  distance  from  the  origin  to  the  point  x(0,h)  on  the  boundary  of 
the  infinite  element.  The  thape  functions  incorporate  the  decay  factor  H  which 
decrease  to  zero  as  5  becomes  Infinite  to  provide  for  finite  Integrals.  The 
factor  H  decreases  at  the  rate  l/r™  as  r  becomes  Infinite. 

At  C  =  0,  the  factor  H  takes  the  value  1  and  tends  to  0  as  5  becomes 
infinite.  These  conditions  are  necessary,  but  not  sufficient,  to  produce  a  good 
approximation.  For  example,  each  of  the  expressions 


1/(1  +  O”  and  r/(R  +  5"*) 


also  satisfies  the  two  preceding  requirements,  but  produces  poor  results  because 
neither  is  a  good  approximation  of  l/r"*  for  moderate  values  of  C.  The  numerical 
Integration  procedure  computes  values  of  the  Integrand  for  moderate  to  large  values 
of  C,  and  the  shape  function  H  must  provide  a  good  approximation  for  these  values. 

Care  must  be  taken  In  forming  the  Infinite  elements  to  avoid  a  mapping  that 
folds  over  on  itself. 3  This  may  happen  If  the  Infinite  elements  are  placed  on 
the  fluid-structure  Interface  and  a  segment  of  the  Interface  Is  not  convex  (as 
shown  In  Figure  2). 


Figure  2  -  Nonunique  Mapping  for  Infinite  Elements 

As  shown  in  Figure  2,  the  shaded  area  will  be  Integrated  over  three  times;  once 
by  each  of  the  Integrations  for  the  center  infinite  element,  and  for  the  elements 
on  either  side  of  the  center  element.  One  method  of  avoiding  this  nonunique  map¬ 
ping  Is  to  make  the  boundary  between  the  finite  element  region  and  the  infinite 
element  region  convex  by  inserting  finite  elements  as  shown  in  Figure  3.  In  this 
case  the  nonunique  mapping  could  also  be  avoided  by  changing  the  direction  of  the 
vectors  nj,  but  the  effect  of  such  a  change  In  the  formulation  on  the  performance 
of  the  infinite  elements  would  have  to  be  determined. 


RESULTS  OF  THE  TEST  PROBLEMS 

The  added  mass  approximation  assumes  that  the  pressure  in  the  acoustic  field 
satisfies  the  potential  equation.  Therefore,  the  pressure  field  in  terms  of  polar 
coordinates  in  the  infinite  exterior  region  is  composed  of  modes  of  the  form 

Pn(r.6)  =  (cos  nO)/r" 

It  would  be  expected  that  an  Infinite  element  would  give  the  best  approximation 
for  a  mode  if  the  rate  of  decay  of  the  mode  matches  that  of  the  element.  To  inves¬ 
tigate  how  well  the  method  approximates  modes  whose  rates  of  decay  do  not  match 
those  of  the  infinite  elements,  natural  frequencies  were  computed  for  modes  n  =  2, 
3,  4,  and  5  using  infinite  elements  with  decay  rates  equal  to  1/r™  for  ra  =  1,  2, 

3,  and  4.  The  frequencies  computed  were  compared  with  frequencies  calculated 
using  an  analytic  solution  for  the  problem.^  Results  from  these  computations  are 
given  in  Figures  4  and  5. 

The  solution  of  the  problem  of  computing  natural  vibration  modes  and  frequen¬ 
cies  for  the  submerged  cylindrical  shell,  modeled  with  finite  and  Infinite  elements 
was  implemented  using  the  NASTRAN  program.  The  actual  model  for  the  problem  repre¬ 
sented  one  quadrant  of  the  plane.  Symmetric  and  antisymmetric  boundary  conditions 
were  used  to  obtain  symmetric  and  antisymmetric  vibration  modes.  To  determine  how 
well  the  solution  had  converged  with  refinement  of  the  finite  element  grid  in  the 
azimuthal  direction,  the  problem  was  solved  using  16  and  32  sectors  per  quadrant. 
For  these  cases,  two  rings  of  finite  elements  were  used  between  the  structural 
elements  and  the  infinite  elements.  The  relative  errors  plotted  in  Figure  4  are 
less  than  3  percent  and  decrease  with  the  refinement  of  the  grid,  although  not 
greatly.  Since  higher  modes  require  greater  detail  in  modeling,  the  errors  become 
larger  for  these  modes.  The  results  show  that  the  solution  is  converged  with 
respect  to  azimuthal  refinement  well  enough  that  futher  refinement  will  not  signif¬ 
icantly  decrease  the  error.  It  can  also  be  seen  that,  although  the  errors  are 
smallest  for  the  modes  whose  rate  of  decay  matches  the  rate  of  decay  of  the  Infinit 
elements,  that  is  if  n  =  m,  the  errors  are  not  extremely  sensitive  to  the  rate  of 


Figure  4  -  Comparison  of  Relative  Errors  for  16  and  32  Sectors 
per  Quadrant,  Using  Two  Rings  of  Fluid  Finite  Elements 


In  the  previous  cases,  two  rings  of  fluid  finite  elements  were  used  between 
the  structural  elements  and  the  infinite  elements.  The  problem  can  also  be  solved 
using  only  infinite  elements  to  represent  the  fluid.  It  is  useful  to  learn  the 
effect  on  the  accuracy  of  the  frequencies  of  the  various  modes  when  only  infinite 
elements  are  used  to  model  the  fluid.  For  these  tests,  16  sectors  of  elements 
were  used  in  the  quadrant.  The  infinite  elements  were  coupled  to  the  structure  at 
the  interface  using  the  same  coupling  technique  that  was  used  with  finite  elements. 
Comparisons  of  solutions  obtained  using  only  infinite  elements  with  solutions 
obtained  using  both  finite  and  infinite  elements  are  shown  in  Figure  5.  The  results 
show  that  there  is  more  sensitivity  to  the  rate  of  decay  when  only  infinite  elements 
are  used.  The  errors  are  small  for  modes  whose  rate  of  decay  matches  the  rate  of 

decay  of  the  infinite  element,  but  become  significant  for  modes  whose  rate  of  decay 

does  not  match  that  of  the  Infinite  element. 

DISCUSSION  OF  RESULTS 

The  results  obtained  in  this  study  Indicate  that  for  fluid-structure  interac¬ 
tion  problems  Involving  structures  submerged  in  large  regions  of  acoustic  fluid, 
one  can  expect  that  the  use  of  infinite  fluid  elements  will  give  good  results. 

Also,  if  results  for  a  broad  band  of  frequencies  are  required,  the  use  of  one  or 
more  layers  of  finite  elements  between  the  structure  and  the  infinite  elements  will 

produce  good  results;  however,  additional  degrees  of  freedom  will  be  introduced 

for  these  elements.  There  is  also  the  possibility  of  using  infinite  elements  that 
have  shape  functions  incorporating  several  rates  of  decay.  These  elements  also 
introduce  additional  degrees  of  freedom.  When  Infinite  elements  are  used  in  one  or 
a  combination  of  these  configurations,  one  can  expect  that  solutions  with  satis¬ 
factory  errors  will  be  obtained. 
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